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ABSTRACT 

Large photometric surveys of transient phenomena, such as Pan-STARRS and LSST, will locate 
thousands to millions of type la supernova candidates per year, a rate prohibitive for acquiring spec- 
troscopy to determine each candidate's type and redshift. In response, we have developed an econom- 
ical approach to identifying SNe la and their redshifts using an uncommon type of optical filter which 
has multiple, discontinuous passbands on a single substrate. Observation of a supernova through a 
specially designed pair of these 'cross-correlation filters' measures the approximate amplitude and 
phase of the cross-correlation between the spectrum and a SN la template, a quantity typically used 
to determine the redshift and type of a high-rcdshift SN la. Simulating the use of these filters, we 
obtain a sample of SNe la which is ~ 98% pure with individual redshifts measured to a z = 0.01 
precision. The advantages of this approach over standard broadband photometric methods are that it 
is insensitive to reddening, independent of the color data used for subsequent distance determinations 
which reduces selection or interpretation bias, and because it makes use of the spectral features its 
reliability is greater. A great advantage over long-slit spectroscopy comes from increased throughput, 
enhanced multiplexing and reduced set-up time resulting in a net gain in speed of up to <~ 30 times. 
This approach is also insensitive to host galaxy contamination. Prototype filters were built and suc- 
cessfully used on Magellan with LDSS-3 to characterize three SNLS candidates. We discuss how these 
filters can provide critical information for the upcoming photometric supernova surveys. 
Subject headings: Supernovae, Dark Energy 



1. INTRODUCTION 

Observations of high redshift Type la Supernova (Riess 
et al. 1998; Perlmutter et al. 1999) reveal evidence of an 
accelerated rate of cosmic expansion due to an apparently 
repulsive-like gravity. Understanding the phenomenon 
of this 'dark energy' provides a fundamental challenge 
to cosmology and new data are vital to this endeavor. 
Initial efforts have focused on measuring the equation of 
state parameter of dark energy, w = P/pc 2 , where P 
is its pressure, and p is its energy density. To determine 
whether the dark energy is a static, cosmological constant 
(w(z) = —1) to higher precision than past supernova sur- 
veys, future ones like the Panoramic Survey Telescope & 
Rapid Response System (Pan-STARRS) and the Large 
Synoptic Survey Telescope (LSST) will attempt to pop- 
ulate the SN la Hubble diagram with 10 4 - 10 6 SNe la 
at redshifts 0.2 < z < 0.8. SNe la remain one of the best 
tools for understanding the properties of dark energy in 
the local volume because they can be discovered in large 
sample sizes and their individual measurement precision 
is high. 

Unfortunately, obtaining a SN spectrum to determine 
its type and redshift is observationally demanding and 
limits the rate at which SNe la can be collected. 1 If Pan- 
STARRS obtains as much time for spectroscopic follow- 

1 The integration time required scales as (l+z) a , with a between 
4 and 6 depending on the redshift and instrumentation. 



up as the SuperNova Legacy Survey (SNLS, Astier et al. 
2006), it will still acquire spectra for only <~ 5% of the 
sample. Purely photometric determination of redshifts 
from multi-band light curves of SNe la (Gong, Cooray 
and Chen 2009), explicitly varying reddening, age and 
light curve shape and assuming the SN is type la, yields 
<r z w 0.1 for 0.2 < z < 0.5 and <r z as high as <~ 1 for 
z > 0.5. The cosmology fitting with only this photo- 
metric data leads to a factor of 4 degradation in the 
dark energy figure of merit compared with the case of 
fitting with spectroscopic data. While the photometric 
approach can be successful in distinguishing core collapse 
SNe from SNe la (Sullivan et al. 2006), the ultimate ac- 
curacy of photometric typing cannot be determined until 
there is a more complete characterization of the varia- 
tions in core collapse SN light curves, specifically how 
often they resemble those of SNe la. An additional com- 
plication arises when using photometric measurements 
to identify SNe la (either for follow-up or for inclusion in 
the Hubble diagram); they cannot subsequently be em- 
ployed for the determination of cosmological parameters 
without introducing a selection bias. Rather, it is impor- 
tant to make use of the wavelengths of spectral features 
(which arc far narrower than broadbands) for type and 
redshift determination, particularly since SN classifca- 
tion is based on spectroscopy (e.g., see Filippenko 1997). 

Here we propose a method of observation that will ex- 
ploit the speed of photometric observations and most of 
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the accuracy of spectroscopic observations to determine 
the types and redshifts of SNe la by using customized op- 
tical filters, each with multiple, discontinuous passbands 
on a single substrate. This approach to characterizing 
SNe la is more economical and can be used to more fully 
harvest the yield expected from large-scale SN surveys. 
In §2, we quantify the fortuitous periodicity in the spec- 
tral features of SNe la that enables one to distinguish 
them from other common transients and determine their 
redshifts when imaged with a semi-periodic passband. 
SN photometry measured with two of these filters ap- 
proximates the cross correlation of a coarse SN la spec- 
trum with a template SN la spectrum commonly used to 
determine the type and redshift of a SN la candidate. In 
§3, we describe our design of these filters and their use. 
We discuss the advantages of these filters as a new tool 
for the follow-up of large-scale supernovae surveys in §4. 
In §5 we present a demonstration of the SNACC filter 
prototype design. 

2. CROSS CORRELATION FORMALISM 

2.1. Spectroscopic Characteristics 

Among common transients (see Fig. 1), SNe la are 
distinguished in the range 3500 — 5500 A by large un- 
dulations in flux which are deep, broad and surprisingly 
sinusoidal. The origin of these spectral features is the 
blending of overlapping P-Cygni profiles of singly ionized, 
intermediate-mass elements (e.g., see Filippenko 1997, 
Baron et al. 1996). SNe la are classified spectroscopi- 
cally by the absence of hydrogen in their spectra and the 
presence of Si II (Ca II, S II, O I, Mg II as well as Fe II are 
also present). Individual spectral lines (absorption and 
emission) have widths of ~ 100 — 150 A, due to expan- 
sion velocities ~ 10000 km s^ 1 . In contrast, the spectra 
of SNe II are continuum-dominated with only a few indi- 
vidual hydrogen P-Cygni profiles because their hydrogen 
shells are intact. The blended spectral features in SN 
la have equivalent widths that are typically 3-4 times 
greater than those in SNe II spectra at 3000 — 5500A 

Although the spectral classification of SNe la was orig- 
inally based on the presence or absence of specific ions, 
the loss of some features outside of the observable wave- 
length range at high redshifts and the typically low 
signal-to-noise ratios (SNRs) of such spectra preclude 
relying on any single ion. As a result, the entire spec- 
trum is considered. Supernova type is determined by a 
strong resemblance to a template spectrum U(X) at a 
trial redshift, z s , where each type is given by a different 
i. Following Blondin and Tonry (2007), one determines 
the most likely template and redshift to be those which 
maximize the cross-correlation function of the template 
and spectrum: 

c{z s ,i) = Js(\)*ti(\x [l + z s ])dX. (1) 

Ideally, high SNR spectra would be obtained for all 
SNe in a survey. In practice, as seen in the ESSENCE 
survey (Matheson et al. 2005) or SNLS survey (How- 
ell et al. 2005, Balland et al. 2009), the need to ob- 
tain many spectra with limited observing time results in 
the collection of many spectra with low SNR degrading 
the information content of the measured cross correla- 
tion function to the point where multiple z s are possible. 



The cross correlation function of both a high and low 
SNR SN la spectrum with SNe la and SNe IIP spectral 
templates is shown in Fig. [2| As expected, the SN la 
templates at the right redshift have a higher correlation 
function amplitude than the SN II templates. For high 
SNR (~ 50 per ~ 10.0 A bin), the correct redshift is 
found with a z = 0.01. However, for spectra with low 
SNR (~ 10 per ~ 10.0 A bin), a periodicity in the SN 
la correlation as a function of trial redshift is apparent, 
favoring a discrete set of preferred redshifts, zx r ue, and 
ZTme ± n x 0.2 (n = ±1, ± 2, . . .), where each local 
peak yields a z = 0.01. With a photometric redshift of 
<j z < 0.08, or a host spectroscopic redshift, this degener- 
acy can be removed. 

We can readily see the cause of this redshift degener- 
acy by viewing the spectra of SN la and SN II spectra 
in Fourier space. S(k) and T(k), the discrete Fourier 
transforms of the supernova and template spectra, can 
be written as, 

JV-l 

S, T(k) =J2 S > t(n)e- 2mk/N , (2) 

n=0 

where we have assumed spectra are observed with N 
equal-sized wavelength bins with s(A) = s(n) and i(A) = 
t(n). n represents the wavelength bin such that A = 
n x kazMii ; anc i = 2ir/n. Following Blondin and 
Tonry, we remove the color information in the spectra 
by dividing it by a cubic spline before taking its Fourier 
transform. The resultant SN la and SN II spectra in 
Fourier space are shown in the top panel of Fig. [3] 
The significant feature of the SN la spectrum is the 
high power at wavenumbers at k — 6, corresponding to 
wavelengths AA ~ 500 A, or spectral feature widths of 
~ 250 A. In the bottom panel of Fig. [5J a sine function 
with this wavelength is compared to a SN la spectrum 
showing their similarity in the range 2500 — 5500 A. This 
periodicity in the spectral features appears to be a for- 
tuitous characteristic of SNe la. While it is possible that 
multiple spectral lines that occur in a narrow wavelength 
range could produce such a sinusoidal shape, randomly 
populating the phase space with a comparable density 
of lines shows that a Fourier transform correlation am- 
plitude as strong or stronger than that seen in Fig. [3] 
occurs < 5% of the time. The distinctive shape of SN 
la spectra is recognized as 'wiggly' by high-redshift SN 
spectroscopists even without certainty of its redshift. 

The correlation amplitude for lower wavenumbers (k < 
5) depends on the intrinsic colors of the SN and redden- 
ing. While such spectral color information is not utilized 
here it may be used elsewhere for distance determina- 
tions. 2 The correlation amplitudes for higher wavenum- 
bers (k > 10) show features that are not detected with 
significance in low SNR spectra. Therefore, from spec- 
tra with low but not atypical SNRs, the only informa- 
tion available to determine the type and redshift of a 
candidate SN la is the power and phase of the k — 6 
mode. The corresponding wavelength scale (500 A) for 

2 In practice, SN colors are more accurately measured from pho- 
tometry from which host light can be subtracted. This also avoids 
chromatic slit losses which occur from long slit spectroscopy ob- 
tained off the parallactic angle (Filippenko 1982). 
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this wavenumber scales with redshift as 500 Ax (1 + z). 
Thus, the peak wavenumber depends on z s such that 
k p = 6/(1 + z s ), but changes little over a modest frac- 
tional change in (1 + z s ), as shown for 0.3 < z < 0.7 in 
Fig. |4j As expected, there is a large correlation ampli- 
tude near k = 4, or 500 Ax (1 + 0.5) w 750 A, throughout 
this redshift range. Thus, a SN la feature centered at a 
true rest-frame wavelength At, and observed at Xo due 
to the redshift zt, can be confused with a similar rest- 
frame feature at Xf such that 

A = A T (1 + z T ) = X F {\ + z T + nAz), (3a) 
where A F - A T = n x 500 A, (3b) 

then 

n x 500 A x (1 + zt) 

Az = — ^ '-. (4) 

At 

For example values of At = 3600 A, zt — 0.5, then 
Az w ±0.2 for n = ±1. Therefore, the undulations in 
SN la flux with periodicity of 500 Ax (1 + z) cause the 
redshift degeneracy of Az ~ 0.2, as in Fig. [2] 

SNe II spectra do not share this periodicity and their 
normalized correlation amplitude at k = 6 is 3.5 times 
lower than that of a SN la, as seen in Fig. [3| Thus, 
measuring the power of SNe la and SNe II spectra at 
k = 6 allows one to distinguish between these two types 
of SNe and the phase of that power provides the redshift 
of a SN la. 

In practice, the determination of SN la type and red- 
shift may be acquired with just two samplings of the 
cross-correlation function at different trial lags. An ex- 
ample of this is shown in the right panel of Fig. [2| 
Two samplings break the degeneracy of the phase, tt, 
of the sine curve and the two measurements are opti- 
mally made with phases 4>\ — 0, (f>2 — 7r/2 (Az as 0.05 in 
redshift space) . Since these two measurements are out of 
phase, at least one of the amplitudes measured reveals 
high power at k = 6 if the SN type is la. If a redshift 
prior (photometric or host spectroscopic) with a z < 0.1 
is also available, one can also confidently determine a SN 
la redshift with a z ss 0.01, similar to the precision of a 
local peak in the SN la cross-correlation function. 

2.2. Imaging Approach 

The convolution theorem states that the convolution 
of an input spectrum and a template is equivalent to the 
product of their Fourier transforms: 

F(s(n) * t{n)) = b x S(k) x T{k). (5) 

Our definition of the cross-correlation, given in Eq. 1, 
is similar to a convolution and thus Eq. 5 is applica- 
ble. Since low SNR spectra have most of their infor- 
mation in the power and phase of the k — 6/(1 + z) 
mode, the template that maximizes the product of the 
fourier transforms, and thus the amplitude of the cross- 
correlation, should have the majority of its power in the 
k = 6/(1 + z) mode as well. Therefore, we can replace 
the SN la template spectra with sinusoidal functions with 
A = 500 x (1 -\-z) with no loss in information in the cross- 
correlation function. 

These simplified templates, t(X,4>i) and t(X,<f>2) 7 can 
now be replaced with optical transmission functions (i.e., 



optical filters) which can be used to evaluate the cross- 
correlation amplitude and phase by direct imaging. We 
now define CC as the SN magnitude through the sinu- 
soidal filter and B as the SN magnitude through a broad- 
band filter. CC — B then measures the cross-correlation 
amplitude as shown in Fig. [2] and can be expressed as 

CC - B-ZP = -2.5 x (log w (^fj) -ZP (6) 

where Fee is the flux through the template filter, Fb 
is the flux through the broadband filter and ZP is a 
zeropoint that will be discussed later. 

Here after we will refer to optical filters which mimic 
t(X,4>) as SuperNovAe Cross Correlation (SNACC) fil- 
ters. With two such filters, these measurements identify 
a SN in the CC\ — B (mag) versus CCi — B (mag) plane. 
If one of these SNACC filters is chosen to be it/ 2 out of 
phase with another filter, then SNe la, over a moder- 
ate redshift range, delineate a circle in this correlation 
space with values of CC — B ~ 0.25 mag when the spec- 
tral features of the SNe are correlated with the template 
bandpasses, and CC — B ~ —0.25 mag when the fea- 
tures are anti-correlated with the template bandpasses. 
The period of the SN la circle with radius 0.25 mag is 
the same as the degeneracy of the cross-correlation ver- 
sus redshift relation for low SNR spectra: Az = 0.2. 
With a host photometric redshift (a z < 0.08), two mea- 
surements with (Jcc-B ~ 0.04 mag in CC — B space 
reveal the redshift of a SN la to a z = 0.01, the precision 
inherent in the SN la correlation curve. 

All SNe II (IIP, ILL, Iln) remain clustered near CC - 
B ~ which provides the means to distinguish them 
from SN la. Since SNe II have low correlation amplitudes 
at k = 6/(1 + z) for < z < 1 in the observed frame 
for a broad range of z, there is little difference between 
the relative flux in a template filter and the broadband 
filter. In addition, the few features of SNe II have lower 
equivalent width with respect to the mean flux than those 
of SNe la, so that if one or two features (at most) align 
with the template filter teeth, the value \CC — B\ is still 
smaller for SNe II than for SNe la. 

3. IMPLEMENTATION 

Optical filters with multiple passbands may be de- 
signed to mimic the spectral features of SNe la and 
used to constrain the type and redshift of SNe. While 
important properties of these filters have already been 
discussed (e.g. separation between the passbands ~ 
500 A x (1 + ZMedian) ) ; the precise determination of the 
filter specifications requires greater rigor. Here the opti- 
mal filters are determined as those which maximize the 
precision of the determination of redshift and SN type. 

We define an extended broadband filter, b' , to transmit 
the sum of flux from Sloan g' and r' so Fy — F g i + F r i. 
For convenience, we define ZP = CCg%v — Bqs,v for Eq. 
6 so that a G8V star has CC — B — ZP = mag. Since 
a G8V star has a thermal spectrum with T = 5310 K, 
and similar colors as the spectra of SNe la at z = 0.5 3 , 
setting ZP = CCgsv — Basv provides a useful reference 
for measurements at {CC\ — B, CCi — B) = (0, 0). Stars 

3 Supernovae have temperatures ~ 10, 000 K . At z = 0.5, then 
10,000 K /(l + z) ~ 6500 K. 
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with the same spectral type but different luminosity class 
will differ by < 0.03 mag. 

The flux Fee through a filter P with discontinuous 
passbands or 'teeth' pi, p 2 . . . is given by 



F, 



CC 



S{X{l + z))P{ Pl ,p 2 ...)d\. 



(7) 



A single passband, or 'tooth', p i: can be described 
with three parameters: a central position (Aj), a width 
(AAj) and a transmission fraction (r^). The set of values 
A, AA,r, for a total of 3 x 2 (filters) x n (passbands), 
can be varied to optimize a merit function and thus the 
design of a filter. 

To construct a useful merit function we first simulate 
a generic SN survey with a few assumptions about the 
properties of SN and their redshift distributions. For 
later optimizations, the survey assumptions are designed 
for a specific survey, e.g., Pan-STARRS-I. Our simulated 
population is composed of 5000 SNe la and 5000 SNe 
II. The SN la spectral templates from Hsiao et al. 2008 
are used for our SNe la spectral template; the spectrum 
of SN IIP 1999em (Baron et al. 2000) is used as a SNe 
II spectral template. SNe Iln and SNe ILL spectra have 
even narrower features than those of SNe IIP, and should 
be even easier to discriminate from SNe la, so SNe IIP 
will represent an upper limit for difficulty in discriminat- 
ing SNe la from SNe II. SNe Ib/Ic/Ibc, as well as other 
transients, are rarer and will be deferred to §3.5. The 
number of simulated SNe grows with the redshift due to 
increasing volume, then decreases with redshift as the 
brightnesses of all SNe no longer exceed the magnitude 
limit. We find the redshift distribution of our simulated 
SN la survey can be approximated as a truncated Gaus- 
sian distribution for a limit of r' = 25 mag and mean SN 
la brightness of —19.16 mag (Richardson et al. 2002): 
z = 0.55, a = 0.15, z mm = 0.30, z max = 0.75. The SNe 
II redshift distribution, with mean brightnesses for SNe 
IIP, SNe IIL, SNe Iln of -16.6, - 17.4, - 18.78 mag, 
can be approximated as z — 0.45, a = 0.15, z m { n — 
0.20, Zmax = 0.7. For each SN, we determine its CC — B 
values from its template spectrum, redshift and the filter 
parameters. We then simulate noisy measurements with 
our assumed photometric errors. 

Because our goal is to find the type of SN, and if a 
la, its redshift, the following chi- squared statistic is min- 
imized with respect to type and z for an individual SN: 

X ^ YPC ,z)J^- CB ^)\(^- CB ^ ' 

V a l J V °2 



+ 



(8) 



where CBi and CB 2 are simulated measurements for 
CC\ — B and CC 2 —B, CB\ Z and CB 2z are the expected 
values of CC — B positions for a SN of a particular type 
and redshift, <j\ i2 are the simulated measurement errors 
of CBi and CB 2 , z p is the available redshift prior and a z 
is its error. The measurement errors, CT12, include image 
calibration errors and photon statistics. We estimate the 
errors resulting from empirically calibrating the CC — B 
value of the supernova from stars in the field to be ~ 0.02 
mag (see §5). The photon statistics are calculated from 
the brightness of the SN and the amount of time alloted 



for imaging on the Low Dispersion Survey Spectrograph 
(LDSS-3, upgraded from LDSS-2, Arlington- Smith et al. 
1994), on the 6.5-m Magellan Clay Telescope. 

We also expect some variation of the SNe la spectra, 
and hence CC-B, which correlate with the luminosity- 
light curve shape parameter or 'stretch', the age of the 
SN and with individual reddening (Guy et al. 2007). 
The remaining, intrinsic variation of SNe la spectra is 
small compared to variations arising from these three 
parameters (Fig. 4 of Guy et al 2007). 

Depending on the circumstances in which the SNACC 
filter data will be analyzed, we may be able to make 
use of knowledge about the age, stretch and reddening 
of a SN la. The available SN light curve may be com- 
plete, sparsely-sampled, or consist of only the discovery 
point. If a well-sampled light curve is obtained, then a 
SN la template with the corresponding age, stretch and 
reddening of the SN candidate would be used to improve 
the chi-squared diagnostic (Eq. 8). Likewise, a z depends 
on whether a host galaxy redshift is known photometri- 
cally (e.g., a z = 0.08), spectroscopically (<r z = 0.01), or 
unavailable (a z = 00). Classification resulting for these 
various follow-up scenarios as well as the effects of age, 
stretch and reddening on the success of determining the 
type and redshift will be discussed in §3.3, §3.4. 

3.1. Merit Metrics for Optimizing Filters 

The two categories of metrics to be used for optimiza- 
tion are the success of determining the type of SN and 
if la, its redshift. Let's assume the two most common 
SN populations, SNe la and SNe II. We define Yj a as the 
fraction of known SNe la that are inconsistent with SNe 
II. Similarly Yjj is the fraction of SNe II that are incon- 
sistent with SNe la. The purity of the SNe la sample 
that would be added to the cosmological sample is then: 



Purity 



Y Ia x N Ia 

Yia x N Ia + (1 - Y n ) x A, 



(9) 



where Nj a and Njj are the numbers of known SNe la 
and known SNe II respectively. 

The success in determining the redshift for SNe la is 
quantified by 

o~z spec —zsNAcc where z S p ec is the known 
redshift and zsnacc is the redshift that minimizes the 
chi-squared diagnostic (Eq. 8). 

3.2. Prototype 3-Tooth Design 

To design prototype filters, we assumed the availability 
of a photometric redshift prior with a z k, 0.08 at 0.3 < 
z < 0.8 and a light curve based age of the SN (which is ±5 
rest-frame days from peak brightness when the CC — B 
mags are obtained). We also assumed integration times 
would be sufficient to reach combined photon statistics 
and calibration errors of ~ a\ 2 < 0.05 mag and that 
the photon statistics would depend on avoidance of the 
skylines in the passbands. 

We first attempted to design a simple filter model with 
just two discontinuous passbands. With filters of this 
kind, the ability to discriminate type of SN is not strong 
over a redshift range of Az — 0.5. Next, we searched for 
a three tooth design for our SNACC filters, which is the 
solution we show here. A filter with four teeth, designed 
by optimizing both the throughput of the filters and the 
discrimination power, will be discussed in §6. 
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A pair of 3-tooth filters, which maximizes the purity 
of SNe la, is shown in Fig. [5] and the expected positions 
of SNe la and SNe II are shown in Fig. |U The widths 
of each tooth is ~ 160 A, but with different transmission 
heights (0.42, 1.0, 0.66), thus these 3-tooth SNACC fil- 
ters have an effective width of ~ 325A, or ~ 33% that of 
a broadband filter. 4 The mean wavelength 



< A >= 



/ XT(X)dX 
J XdX 



(10) 



of the two SNACC filters (< A >= 5750, 5650 A re- 
spectively) and the broadband filter (< A >= 5550 A) 
are quite similar which greatly limits the effect of red- 
dening on the observed CC — B values. Also, the band- 
passes of the second filter avoid all major skylines, but 
because the center bandpass of the first filter is near the 
5887 A Na D line, this filter transmits 20% more sky 
flux than the first. As shown in Fig. [fjj SNe la CC — B 
mags trace roughly concentric circles tracking counter- 
clockwise with increasing redshift. As expected, SNe II 
lie closer to the origin due to their lack of spectral fea- 
tures as strong as SNe la, with an average separation 
from the origin of ±0.05 mag. For measurement errors 
of cr m = 0.04 mag (as we obtained with the first use of the 
SNACC filters, discussed in §5.2), the SNe la and SNe 
II are separated in CC — B space by ~ 3.0 — 5.0 x a m . 
The actual fraction of SNe la correctly classified, Y] a , is 
found to be 99.0% and the purity of SNe la, assuming 
equal sizes of SNe la and SNe II populations, is 98.0%. 
Measurement errors of 0.03 and 0.05 mag, which we con- 
sider lower and upper limits for our measurements, result 
in purities of ~ 99.0% and ~ 95.0% respectively. 

SNe ILL (using spectral templates that are updated 
versions of those given in Nugent et al. 2002 5 ), like SNe 
IIP, have an average distance from the origin of CC — B 
space of « 0.05 mag. SNe Iln spectra are mostly con- 
tinuum, with few narrow lines that form in low velocity 
circumstellar matter. Since the width of these lines is so 
small, the average distance of SNe Iln in CC — B space 
from the origin is f=a 0.01 mag. 

The histogram of Az = z spec — z snacc for all the SNe 
la is shown in Fig. [7] assuming a host photometric error 
er 2 = 0.08 and measurement errors — 0.04 (mag). 
The redshift distributions are centered around Az = 
with a standard deviation around the center of a z = 0.01. 
Due to the cross-correlation redshift degeneracy and the 
limited assumed precision of the host photometric red- 
shift, there are over-populated tails at Az — ±0.2 that 
hold « 3% of the distribution. Without a photomet- 
ric redshift, the tail population is rs 45%. An error in 
z = ±0.2 should be easily distinguishable in the Hubble 
diagram since the difference in distance for two SNe la 
with this redshift separation is ~ 1.0 mag. 6 

3.3. Effects of Light Curve Shape, Age and Reddening 
on CC — B magnitudes 



4 Elsewhere (Scolnic et al. 2010) we present a 4-tooth design 
which a includes a 4th tooth between 1st and 2nd as seen in Fig. 
[H] and achieves a width of 600 A without loss in purity. 

5 http : / /supernova, lbl . gov/~nugent/nugent_templates .html 

6 Ad = 51og^ _ 51og c(zr + a2) = 51og^L « 



1.0 mag for zt ~ 0.5. 
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Variations in SN la spectra are mainly due to varia- 
tions in reddening, light curve shape and age. The ef- 
fects of ignorance of these parameters (using the SNe la 
variations modeled in Guy et al. 2007) on our purity and 
redshift errors are shown in Fig. [8] and Fig. [9] respec- 
tively. Typical usage of SNACC filters would avoid these 
errors using light curves to choose a SN la spectral tem- 
plate corresponding to the appropriate age, light curve 
shape and reddening. 

Due to the near-equality of the mean wavelength of the 
SNACC filters and the broadband filters, the fraction 
of SN la correctly classified at most redshifts is quite 
insensitive to reddening, affecting Yj a at the 1% level, as 
shown in Fig. [8] Selection purity is also insensitive to the 
light curve - luminosity relation, varying by ~ 2% (except 
for z > 0.75). Ignorance of the SN age has a stronger 
impact on Yj a and limits the time range in which it is 
optimal to use SNACC filters for follow-up. The purity 
is very robust for ages of —7 to ±5 days but degrades for 
older SNe la. Y Ia is < 90% for SNe with z < 0.35 past +5 
days and for SNe with z > 0.75 past ±8 days. For SNe la 
with 0.35 < z < 0.75, Yj a on average is ~ 95% up to an 
age of 10 days, and even for ages past this, Yj a for parts of 
the redshift range remains as high. Prior to 5 days before 
peak, the overall purity remains > 95%, but SNe la with 
z ~ 0.55 (for —10 days) or z ~ 0.65 (for —15 days) will 
approach the origin in CC — B space. While it is difficult 
to simulate how peculiar SNe la with similar early ages 
and high redshifts appear in CC — B space because of 
the limited set of UV spectra, SNe like SN 199 IT (Jeffrey 
et al. 1992), for redshifts 0.25 < z < 0.40, are located 
in CC — B space with the same phase and at most a 
~ 0.05 mag relative distance to the origin as their normal 
SNe la counterparts (and Yi a /a z values of 92%/0.012). 
For these peculiar SNe la, the strengths of the spectral 
features vary from those of normal SNe la, but since three 
separate passbands are correlated with the spectrum, the 
overall difference between these two categories of SNe la 
in CC — B space is not large. The effects of age, light 
curve shape, and reddening are most noticeable at the 
limits of our redshift range (z ~ 0.35,0.75) since the 
spectral features at these redshifts that correlate with the 
passbands are weak compared to the correlating features 
at z ~ 0.5. 

Both reddening (0.0 < E(B - V) < 0.4) and stretch 
(0.85 < s < 1.15) cause errors in \zsnacc — zs pe c\ that 
are all < 0.02 as shown in Fig. [9] . Anytime before 5 
days past peak, age has a limiteaeffect on the redshift 
errors (<~ 0.02), but afterwards the effect is much more 
significant (typical errors ~ 0.05) for all redshifts. 

3.4. Usage Scenarios 

SNACC filter measurements may be used in conjunc- 
tion with other data to yield more precise results. For 
example, the host galaxy of the SN may be apparent, and 
the multi-band light curves may have been acquired. If a 
host galaxy is visible, a photometric redshift (a z — 0.08) 
or spectroscopic redshift (a z — 0.01) may be obtained. If 
a light curve has been acquired, values of the age, stretch, 
reddening and photometric redshift (cr z /(l + z) = 0.06, 
based on Sullivan et al. 2006) can be fit on the con- 
dition the SN is la. An available photometric or spec- 
troscopic redshift breaks the Az = 0.2 cross-correlation 
degeneracy of SNe la and knowledge of the age, stretch 
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and reddening reduces the spectral variations because 
the appropriate spectral template can be used. 

As realistic examples, the Medium Deep Pan-STARRS 
survey, which images the same area of sky every four 
days, should provide well-sampled light curves. In con- 
trast, the 3tt Pan-STARRS survey covers the entire sky 
only twice a month, and will not provide light curves 
well-sampled enough to fit accurate values of the type, 
stretch, reddening or age of the SN. 

The combinations of possible redshift priors (photo- 
metric, spectroscopic, none) and light curve based priors 
(yes, no) make up six distinct scenarios. Table I presents 
these scenarios and the results expected after SNACC 
follow-up. We simulate a realistic set of SNe la with a 
uniform distribution of ages within ±5 rest frame days 
from peak, a gaussian distribution of stretches with s = 1 
and a s — 0.1 and an exponential distribution of redden- 
ing values, with a peak at E(B — V) = and r = 0.1. 

Table I shows that the fraction of SNe la correctly clas- 
sified and the redshift determination both improve with 
the amount of information one obtains about the super- 
nova or host galaxy. With a well-sampled light curve 
and host galaxy but no SNACC filter measurements, we 
assume a baseline of Yj a — 80% and a z /(l + z) w 0.06. 
Yj a is 97.0% if one obtains a well sampled light curve of 
the SN because then one can use the correct template 
in the chi-squared diagnostic of Eq. 8. Similarly, the 
precision of the redshift determination is also greater if 
a well sampled light curve is obtained. The best case is 
the scenario in which one obtains a spectroscopic redshift 
from the host galaxy in addition to a well sampled light 
curve of the SN. Here, Y] a is 99.0%, although this is not 
significantly better than the case in which one obtains 
only a photometric redshift of the host. Generally, if one 
gathers supplementary information about the SN in ad- 
dition to the SNACC filter measurements, one achieves 
Y Ia >~ 95% with a z < 0.015. 

3.5. Other Transient Candidates 

Until now, only the properties of SNe la and SNe II 
in CC — B space have been considered because they are 
the most common and therefore their separation is most 
important to good SN la selection. There are other types 
of transients, and we list here some that may be confused 
with SNe la (Matheson et al., 2005) and their character- 
istics in CC — B space. 

SNe Ib/Ic/Ibc The progenitors of SNe Ib/Ic/Ibc 
(hereafter called SNe Ib-c) are believed to be stars 
stripped of their original H-rich envelope (Georgy 
et al. 2009) and the spectral features of these su- 
pernova are similar to those of SNe la (as seen in 
Fig. 1). The SNe Ib-c spectral templates we use 
are updated versions of those given in Nugent et 
al. 2002. The characteristic spectroscopic differ- 
ence between Type la and Ib-c SNe is the deep 
absorption trough at 6150 A in Type la spectra, 
which is due to the blue-shifted Si II AA6347, 6371 
feature (Homeier 2005). Unfortunately, this is red- 
shifted out of the optical passband for z > 0.4, 
making it difficult to distinguish SNe la from SNe 
Ib-c. Although there are other spectroscopic dif- 
ferences between SNe la and Ib-c (e.g. Si II A4130, 
Coil et al. 2000), they are too slight for the two 



types with the same redshift to appear in different 
locations in CC — B space. There are two branches 
of redshifts, 0.1 < z < 0.18 and 0.30 < z < 0.38 in 
which the SNe Ib-c are clustered in the center, but 
for other redshifts, their CC — B positions are often 
less than 0.05 mag away from a SN la with the same 
redshift. However, they are rare enough (< 3% 
of photometric SN la candidates targeted for spec- 
troscopy in the ESSENCE project (Matheson et al. 
2005)) to not simultaneously affect our SN la pu- 
rity. We also note that SNACC filter observations 
of SNe Ib-c fair no worse than typical high-redshift 
spectroscopy whose SNR is frequently (Riess et al. 
1998, Matheson et al. 2005) too low' to identify Si 
A4130 and which is not even present for many SNe 
la. For the set of ESO/VLT spectroscopy of SNLS 
candidates (Balland et al. 2009), 38 out of 124 of 
the classified SN la spectra were called 'probable' 
SNe la because other types, in particular SNe Ic, 
could not be excluded given the SNR or the phase 
of the spectrum. In any case, one learns the correct 
redshift for SNe Ib-c which aids their elimination 
from the Hubble diagram with Bayesian techniques 
(Kunz et al. 2006). 

AGN There are several subclasses of AGN (e.g. Seyfert 
1, Seyfert 2, LINER) and the AGN in each subclass 
have unique spectral features (Osterbrock 1989). 
Here we use the set of 12 AGN spectra obtained 
by the ESSENCE project which are particularly 
relevant because they were photometrically iden- 
tified as SN la candidates. In the project, AGN 
were discovered with redshifts 0.18 < z < 2.6, so 
we must include a large redshift range when deter- 
mining the location of AGN in CC — B space. For 
the ESSENCE set of spectra, over the redshifts in 
which the spectra cover the entire optical range, the 
AGN will deviate from the origin > 0.05 mag only 
~ 5% of the time. Thus, > 11 out of 12 ESSENCE 
AGN would have been rejected from the SN la sam- 
ple. For spectra at select redshifts (e.g. spectra 
similar to those of Seyfert 1 AGN at z <~ 2.57, see 
Francis et al. 1991 for spectral templates), the po- 
sition of AGN in CC — B space may deviate as far 
from the origin as SNe la. Many AGN can be dis- 
criminated from SNe la with broadband color mea- 
surements and possibly a time-series history, and 
SNACC filters would provide an additional method 
to differentiate the AGN from SNe la. 

Variable Stars During the course of the ESSENCE 
project, as well as other recent supernova surveys, 
several variable stars as well as unclassified stars 
were discovered. In the spectra published from the 
ESSENCE project, the spectra of the stars were 
mostly continuum and are located near the origin 
of CC — B space with little deviation. Thus, rarely 
should any variable stars be confused with SNe la. 

4. OBSERVATIONS 

4.1. SNACC Filter Fabrication for LDSS-3 

The technology used to fabricate SNACC filters, with 
multiple discontinuous passbands on a single substrate, 
has become popular for a wide range of applications in 
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many fields, like fluorescence microscopy in biology (e.g. 
Yelin and Silberberg, 1999). We canvassed multiple ven- 
dors but ultimately Iridian Spectral Technologies Ltd. 7 
was best able to match our specifications for the proto- 
type within a reasonable cost. 

Iridian constructs Thin Film Filters (TFF, see Macleod 
2001, for a review) by layering materials of alternating 
refractive index with thicknesses on the order of visual 
wavelengths, ~ 5000 A, to produce optical filters that 
transmit and/or reflect specific wavelengths. A thin- 
film coating is a stack of such layers having boundaries 
of greater refractive index (reflecting light 180° out of 
phase) and less refraction index (reflecting light in phase) 
that if at precise spacings will modify the reflected and 
transmitted components by interference. 

The filter design team at Iridian use a proprietary 
design and process control technology, PrecisionSpec- 
trum@IST, to develop the optimal multilayer design for 
the requested filter specifications and requirements (both 
physical and optical). The SNACC filters were designed 
with a TFF layering on the front-side and simple two- 
layer anti-reflective type coating on the backside. 

The filters were specified for use in the Low Dispersion 
Survey Spectrograph (LDSS-3, upgraded from LDSS-2, 
Allington-Smith et al. 1994), on the 6.5-m Magellan Clay 
Telescope. Our specified transmission function of the fil- 
ters was adjusted to include the effect of the camera ef- 
ficiency. Our request, Iridian's expected des ign, and the 
design they achieved are shown in Fig. |10[ The band- 



passes of the filters Iridian produced had width, height 
and position values that were all within 5% of our spec- 
ifications. 

In the LDSS-3 setup, an input f/11 beam enters a colli- 
mator before passing through a SNACC filter. The light 
is then focused by the camera onto an external detector 
with a final focal ratio of f/2.5. Due to off-axis rays, 
there is a correlation between the incident angle on the 
filter and the radial field position. For 2/3 of the field, 
there is a maximum angle of incidence of ~ 4.6° of the 
beam on the filter. Because of the construction of the 
interference filter, there is a shift in the effective filter 
wavelength, to the blue (i.e. at 6000 A, this is a shift of 
about 10A). The filters were analyzed in the optics lab 
at Johns Hopkins University to measure their transmis- 
sion functions in a similar setup with a tungsten halogen 
lamp, silicon photodiode, and collimator. We measured 
filter transmissions which were very close to what Irid- 
ian stated they achieved in a parallel beam, and that the 
out of band rejection was > 99.5% over the wavelength 
3000 — 11000 A. When collimated beams have an angle 
of incidence of up to 5° on the filter, the shifts in the 
bandpasses were fairly coherent with values up to 8 A, 
negligible for our purpose. 

4.2. SNACC Filters Imaging and Spectroscopy with 
LDSS-3 

To test the utility of the filters for SN classification, we 
observed three supernovae and spectrophotometric field 
calibration stars. Several supernova candidates were dis- 
covered in the SuperNovae Legacy Survey (SNLS, Astier 
ct al. 2006) with the Megaprime camera at the Canada- 

7 For more information see http://www.iridian.ca/ 



France-Hawaii telescope after observations on 2008 May 
1 (previous observation date was 2008 April 14) in the 
SNLS D2 field (RA 10:00:29, DEC 02:12:21, 1 square de- 
gree). Along with the candidates, SNLS provided the 
candidates' photometric redshift if it was assumed they 
were SNe la using the SN's broadband color measure- 
ments, g,r,i,z measurements around peak, and photome- 
try of the host. A significant advantage of the SNACC 
filters over long-slit spectroscopy is that they easily al- 
low for host light subtraction. In this case, however, we 
chose to follow-up reasonably isolated SNe since we did 
not have time to obtain SNACC filter host templates. 
We chose the candidates 08D2nu (NU) and 08D2oe (OE) 
because their available (though minimally sampled) light 
curves exhibited the pre-peak rise and appropriate g — r 
colors of SNe la light curves (sensitive to reddening) , and 
08D2jh because its light curve exhibited the post-peak 
plateau of a SN IIP light curve. Referring back to §3.4, 
the follow-up scenario in which our observations of these 
candidates took place is the one in which a light-curve 
is sparsely populated and there is a photometric red- 
shift from the host galaxy. The calibration stars, a G8V 
star (r' = 16.7 mag) and a M4V star (r' = 17.4 mag), 
were selected from the Sloan Digital Sky Survey Data Re- 
lease 7 (SDSS-DR7) (Abazajian et al., 2009), which pro- 
vided photometry and spectra. The stars' types were de- 
termined from cross-correlating their spectra with spec- 
tra from the Gunn & Stryker spectrophotometric atlases 
(Gunn & Stryker, 1983). The G8V star was chosen be- 
cause its CC — B values should be near the origin, and 
the M4V star, with its strong spectral features, was cho- 
sen because it should deviate as far from the origin as 
SNe la in CC — B space and would be helpful to verify 
the transmission function in the beam. As we will discuss 
later, the CC — B values of the G star (with g — r < 1) 
only depend on its luminosity class because of the rela- 
tion between luminosity class and g — r color, whereas 
the CC — B values of the M star (with g — r > 1) de- 
pend strongly on its luminosity class because of the star's 
strong spectral features. 

The observations obtained for the SNe and calibration 
stars are given in Table 2 and 3 respectively. We observed 
SN NU and SN OE on UT 2008 May 3 with the SNACC 
filters on LDSS-3 for 300 seconds each, and SN JH on 
May 5 for 300 seconds as well. The G8V and M4V stars 
were observed for exposure times of 180 seconds on UT 
2008 May 3. Flat-field images and bias frames were also 
acquired. A custom pipeline was used to perform pho- 
tometry of the bias-subtracted, flat-fielded images, using 
aperture photometry via the APER routine in IDL and 
point-spread function photometry via the DAOPHOT 
routines in IDL. 

We also acquired spectra of the SNe NU, JH, OE on 
UT 2008 May 4, 5, 9 with the VPH-A11 grism on LDSS-3 
with a 1" slit for integration times of 2 x 1200 s, 2 x 900 
s, and 2 x 1800 s respectively. The VPH-A11 grism (400 
lines/mm) allowed us to acquire 4000 - 9900 A spectra 
with an average resolution A/ A A w 860 across the en- 
tire chip. Dispersion along the chip was 1.9 A /pixel. 
All spectral observations were accompanied by HeNeAr 
arc lamps exposures to measure dispersion and the pixel 
to pixel response was removed with flats from a quartz 
lamp. The spectral data were reduced using the standard 



8 



IRAF 8 packages including APALL (extraction), DISP- 
COR (assigning the dispersion function from the cali- 
bration arc lamp) and CALIBRATE (using sensitivity 
function to flux calibrate the spectra). Flux calibration 
(using STANDARD) of the spectra was performed by 
means of spectrophotometric standard stars observed at 
similar airmass on the same night as the SN. For all the 
spectra observed, the slit was always aligned along the 
parallactic angle to avoid differential chromatic refrac- 
tion. 

5. ANALYSIS OF OBSERVATIONS 

5.1. SNACC Filter Calibration: Photometry for C and 

M' stars 

Using spectra from the Gunn & Stryker spectropho- 
tometric atlases (Gunn & Stryker, 1983) and the known 
transmissions of the broadband and SNACC filters, we 
use the relation between the CC — B values and g — r 
values for standard stars to determine their zeropoints. 

The CC — B values of the stars are linearly rela ted to 
their g — r values, as seen by the diamonds in Fig. 11|12 
such that one can write the relation for stars 



CC - B = m(g-r) + Y for g - r < 1.0 



(11) 



where CC — B has already been offset by the zeropoint as 
defined previously, ZP = CCqsv — -Bgsv ■ The slope and 
intercepts can be found by a least-squares fit: (mi, Yi) = 
(-0.14, 0.04), (m 2 ,Y 2 ) = (-0.19,0.07). These relations 
are highly linear for g — r < 1.0 and have a dispersion 
of less than 0.02 mag. Stars with g — r > 1.0 generally 
are M-type stars with strong spectral features that break 
the CC — B/g — r relation. 

For a set of n observed stars, each with CC, B, g and r 
values (5 and r have been corrected for zeropoint offsets), 
the CC — B color zeropoints ZPi and ZP2 for a field can 
be determined from: 

CCi - B - ZPi = -0.14( 5 -r-) +0.04 (12) 
CC 2 -B - ZP 2 = -0.19(5 -r) +0.07 (13) 

To find the zeropoints, one minimizes \ 2 with respect 
to ZP 

" (ZP-((CC 4 -i? l )-m(.g 4 -r t ) + r) 2 n 
X = > . j (I 4 ) 



i=0 
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where <?i — < 1.0 and cr, is the quadrature sum of star 
i's measurement errors. Having determined ZPi^, the 
CC — B values of the SNe in the field can be placed in 
our diagnostic CC — B space. 
The results of this calibration method for the G8V field 



and M4V field are shown in Fig. [TTJand Fig. [12} The ZP 
values found for the G8V field were ZP X = 0.33+0.01 and 
ZP 2 = 0.24 ± 0.005. The ZP values found for the M4II 
field were ZP X = 0.33 ± 0.015 and ZP 2 = 0.22 ± 0.01. 

The CC\ — B and CC 2 — B values of all the stars 
in the G8V field and M4V can be combined, as can be 
seen in Fig. 13 and Fig. 14 and we can compare the 



values of the TT8V and M4V stars with their synthetic 

8 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



spectrophotometric measurements (spectra from SDSS 
DR7). To determine the total photometric error, the 
calibration error of the stars were added in quadrature 
with the photon statistics error. The G8V and M4V 
stars were ~ 17 mag in the r' band, and the total photon 
statistics error from the SNACC filters and broadband 
filters was ~ 0.015 mag. As seen in Fig. [13] and Fig. 
14 the synthetic measurements are within one standard 
cTeviation of the photometric measurements. 

5.2. Results of Imaging and Spectroscopy 

The extracted SN spectra, as observed with the LDSS- 
3 spectrograph, are presented in Fig. |15[ Bluer than 
4300 A, the SNR< 5 due to the decreased sensitivity 
of the spectrograph and lower object flux with respect 
to the sky. Template spectra of the same type and age 
were used to graft this part of the spectrum. The SNID 
program (Blondin and Tonry, 2007) was used to deter- 
mine that NU and OE are both SNe la, with redshifts 
z = 0.56+0.01 and z = 0.42+0.01 respectively, and JH is 
a SN IIP with z = 0.23 + 0.01. The quality of correlation 
is determined by the rlap quality parameter, which is the 
product of the correlation height-noise and the spectrum 
overlap parameter, and is 0.63, 0.68 and 0.74 for NU, 
OE and JH respectively. Using the extracted spectra 
as well as transmission functions of the SNACC filters 
and broadband filters, the synthetic spectrophotometric 
CC—B values of the SNe can be evaluated and compared 
to the observed values. 

Details of the SNe ages, broadband colors, and red- 
shifts (host galaxy photo-z, SN photo-z from broadband 
colors, and spectroscopic) are given in Table 4. Obser- 
vations of the SNe were sky-noise limited, so the pho- 
ton statistics error is determined from the dispersion of 
the photometry when a model psf with the source's am- 
plitude is placed at random locations in the field. The 
300 seconds of imaging yielded SNRs of about 35 for 
NU and JH, and 30 for OE. Using the calibration ap- 
proach discussed in §5.1, with Megaprime g' and r' fil- 
ters and SNLS calibrated stars to calibrate the G and R 
values, we obtained calibration errors for (ZPi,ZP2) of 
(0.01,0.015), (0.01,0.14), (0.007,0.009) for NU, OE and 
JH respectively. The positions of the SNe candidates in 
CC — B space are shown in Fig. [16] The photometric 
positions in CC — B space are within one standard devi- 
ation from the corresponding synthetic positions, giving 
us confidence we are properly quantifying the SN flux 
through the SNACC filters with our photometry. 

To determine the probability function a SN is a certain 
type and redshift, we employ the x 2 diagnostic of Eq. 8. 
We make use of the sparse SN light curves to determine 
that our near-maximum SN templates are appropriate 
for the \ 2 diagnostic. The photometric redshifts from the 
host galaxy are 0.54+0.13 for NU, 0.47+0.20 for OE and 
0.27 + 0.03 for JH. The probability distributions for each 
SN are shown in Fig. [17] For SN NU, the probability 
that the SN is type la with z = 0.56 is ~ 20 x the proba- 
bility of any SNe II at any redshift. There are negligible 
probability tails at z = ±0.2 because the photometric 
redshift from the host galaxy is a z = 0.13 and because of 
the NU's CC — B position. The measured redshift error 
uncertainty is thus cr z = 0.01 and the most probable red- 
shift is within this error from the spectroscopic redshift. 
For SN OE, the SN fell in the region of CC - B space 
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where there is a local degeneracy of CC — B values. The 
measurement strongly indicates that SN OE is type la, 
but with a poor photometric redshift z — 0.47 ± 0.2, it 
is difficult to determine the true redshift. The spectro- 
scopic redshift 0.42 can be seen in one of the degenerate 
probability tails. As seen in Fig. [16] the photometric 
and synthetic position are within one standard deviation 
of each other, so a more accurate measurement would 
not improve the classification. While our 3-tooth pro- 
totype suffered from the degeneracy in the upper right 
quadrant of the CC — B space, we were able to correct 
this in our 4-tooth solution, which will be used in the 
future and has no degenerate areas. Lastly, for SN JH, 
the probability that the SN is a SN la at a given redshift 
is < 0.01. Although one can not say the redshift of the 
SN, the low probability for all SNe la would allow us to 
discriminate this SN from SNe la and remove it. 

6. DISCUSSION 

While the use of multi-bandpass filters to measure a 
complex cross-correlation function is novel, such filters 
have had a prior use in astronomy. They were used to 
monitor Blazars in multiple colors simultaneously (Wu 
et al. 2007). The light passing through the different 
passbands may be differentially refracted by an objec- 
tive prism and can then be spatially separated so there 
are multiple color images for each object obtained si- 
multaneously. This system enables accurate tracing of 
contemporaneous color change of the blazars. Surpris- 
ingly, because of the shape of certain Blazar spectra, the 
location of the bandpasses are very similar to those of 
the SNACC filters. 

For our prototype 3-tooth design, we did not attempt 
to minimize the exposure time by maximizing the filter 
throughput. We have begun work on this optimization 
and find that we can create 4-tooth filters, as shown in 
Fig. 18 with effective widths of around 600 A, which is 
than half the width of broadband filters. As seen in 
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Fig 



19 this design avoids any local degeneracies, while 



still maintaining a purity of SNe la greater than 98% 
and a redshift error <j z = 0.01, equivalent to what we 
can achieve with the prototype design. 

Because of gains in throughput, multiplexing and set- 
up, using SNACC filters is greater than an order of mag- 
nitude faster than long-slit spectroscopy. For example, 
with SNACC filters designed for SuprimeCam at Subaru 
(Miyazaki et al. 2002), acquiring images of SNe (SNR 
> 30) is ~ 4 times faster than acquiring spectra (with 
SNR ~ 30, using the Subaru Exposure Time Calcula- 
tor) , and the expected space density of SNe la in a Pan- 
STARRS-like survey allows one to collect 4 to 8 SNe per 
pointing (with a field of view of 0.25 square degrees), 
so using SNACC filters is up to ~ 30 times faster than 
spectroscopy. 

Spectra of high-z SNe (z > 0.5) with a sufficient SNR 
for identification may be acquired on a large telescope 
(e.g., the Very Large Telescope (VLT), 8 m, see Appen- 
zcllcr et al. 1998) in < 1 hour, but the field of views for 
these telescopes' spectrographs is too small (e.g., FORS1 
on VLT, 36 square arcmin) to allow any multiplex ad- 
vantage from multi-slit spectroscopy. On smaller tele- 
scopes (e.g., Magellan Baade and Clay Telescope 6.5 m, 
see Dressier 2004), it is more common to obtain spectra 
of SNe with redshifts 0.2 < z < 0.5 (e.g., Foley 2009), 



and while the field of view may be larger (e.g., IMACS 
on Baade, 239 square arcmin), the viewing depth is not 
large enough to allow any significant multiplex advantage 
from multi-slit spectroscopy. 

In future surveys, measurements with SNACC filters 
should yield a likelihood distribution for each individual 
SN candidate being type la within a measured redshift 
range. Such a distribution is well incorporated into a 
Bayesian cosmology analysis. While there are small tails 
in this probability distribution from the SNACC mea- 
surements of Az = ±0.2, an error in z = ±0.2 should be 
easily distinguishable in a Bayesian cosmological analysis 
of the Hubble diagram since the difference in distance for 
two SNe la with this redshift separation is ~ 1.0 mag. 
Also, while there are certain transients, like SNe Ib-c, 
that are easily mistaken for SNe la, measurements with 
the SNACC filters will still provide the correct redshift of 
the SNe Ib-c. With the right redshift, the rarity of these 
objects, and the known brightness difference between su- 
pernovae types, a Bayesian cosmological analysis should 
easily be able to cope with these contaminants (Kunz et 
al. 2006). 

Aside from speed, there are a number of advantages 
provided by the SNACC filters, which can be valuable 
depending on the application. The measurements with 
SNACC filters are much less sensitive to reddening and 
allow for a more precise determination of redshifts than 
usual photometric methods. Image subtraction allows for 
a measurement free from the host contamination afflict- 
ing spectroscopy, which can be essential near the centers 
of hosts. Indeed, a well known selection effect pertains 
at high redshift by the inability of long-slit spectroscopy 
to identify SNe la for cosmological analysis. This bias 
may be lifted with SNACC filters. And lastly, using the 
SNACC filters is independent of the color data used for 
subsequent distance determinations so there is no selec- 
tion or interpretation bias (Kessler et al. 2009), and 
because it makes use of the spectral features its reliabil- 
ity is greater. Although there are other approaches, like 
the use of multiple Dichroic stacks, that may be able to 
achieve the same goal of fast spectroscopy, SNACC fil- 
ters are far more accessible and economical because they 
can be deployed in available filter slots. 

7. SUMMARY 

To obtain the needed type and redshift information 
of candidate SN, we have developed a method of super- 
nova observation which exploits the speed of photometric 
observations and much of the accuracy and precision of 
spectroscopic observations to determine the type of SN, 
and if la, its redshift. Our approach is based on the use 
of multiple narrow passbands in a single filter. Obser- 
vations of supernovae through two of these filters effec- 
tively provide equivalent information as the evaluation 
of a spectral cross-correlation, typically used to measure 
the redshift and type of a SN from a low SNR spectrum. 
Simulating the use of SuperNovae Cross Correlation Fil- 
ters (SNACC) filters as a follow-up tool, we can obtain 
a sample of SNe la which is ~ 98% pure with redshifts 
of individual precision <r z = 0.01 at a rate up to ~ 30 
times faster than typical high-z SN spectroscopy, a speed 
advantage which could greatly benefit future wide-area 
transient surveys. 

Prototype SNACC filters were built and tested on 
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Magellan with LDSS-3 and used to follow-up three SNLS 
candidates. We acquired images with the SNACC filters 
as well as spectra of each candidate and we were success- 
ful in classifying the three candidates. In the future we 
plan to apply this new technology for the next generation 
of supernova surveys and try to improve the precision of 
the supernova sample and its constraining power for dark 
energy. 
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TABLE 1 

Usage Scenarios 



Scenario 


Details 


Results 








SNe la Correctly Classifed (%) 





Host Galaxy, Light Curve 
SNACC, Host Galaxy w/ Photo-z, no LC 
SNACC, Host Galaxy w/ Spec, z, no LC 
SNACC, No Host Galaxy, w/ LC 

SNACC, Host Galaxy w/ Photo-z, w/ LC 

SNACC, Host Galaxy w/ Spcc-z, w/ LC 



cr z =0.10, <r z /(l + z) a 0.06 a , 

cr z =0.10 

cr z =0.01 

cr z /(l + Si 0.06 a , 
known age, stretch, reddening 
cr z =0.10, <7 Z /(1 + z) Si 0.06 a , 
known age, stretch, reddening 

<7 Z =0.01, 

known age, stretch, reddening 



80.0 (5.00) 
94.5(0.2) 
95.6 (0.3) 
97.0 (0.2) 

98.2 (0.1) 

99.0 (0.1) 



0.090 (0.010) 
0.015 (0.003) 

0.014 (0.002) 

0.012 (0.002) 



These photometric rcdshifts, derived from the broadband colors of the SN, are dependent on whether the SN is known to be a SN 



la. 



TABLE 2 
SN Observation Log. 



Source 



Coordinates [J2000.0] of SN Imag./Spec. UT Date 



(s) Offset [arcsec]" Airmass 



SN 08D2nu 09 h 58 m 32 3 .778 + 02°36'05".98 
SN 08D2oe 10 ,1 01 m 14M36 + 02°04'16".99 
SN 08D2jh 09 h 59 m 31 s .661 + 02°14'51".93 



Imag. 
Spec. 
Imag. 
Spec. 
Imag. 
Spec. 



2008-05-03 300 0.54E 2.03N 1.214 

2008-05-04 2 x 1200 1.220 

2008-05-03 300 ... b 1.434 

2008-05-09 2 x 1800 1.170 

2008-05-05 300 0.01W 0.18N 1.170 

2008-05-04 2 x 900 1.460 



a SN positional offset from nucleus of the host galaxy 
b SNR< 5 of host galaxy for all data 



TABLE 3 

Spectrophotometric Standards Imaging Observation Log. 



Source" 


Spectral Type 


UT Date 


tint 


Airmass 


g' 


r' 








(s) 




(mag) 


(mag) 


SDSS J141101. 76+010418.4 


G8V Star 


2008-05-03 


180 


1.490 


17.19 ±0.01 


16.67 ±0.01 


SDSS J125857.35+145019.7 


M4V Star 


2008-05-03 


180 


1.157 


18.97 + 0.01 


17.40 ±0.01 



The coordinates arc given inside the source name, in HHMMSS of J2000. 



TABLE 4 
Main SNe Characteristics 



SN 


Type 


Age° 


g' 


r' 


Host Gal. Photo-z 6 


SNLS photo-z 


Spec. z d 






(days) 


(mag) 


(mag) 




(if SN Ia) c 




NU 


la 


+1 


23.4 


22.5 


0.54 + 0.13 


0.42 


0.56 


OE 


la 


-3 


23.3 


22.9 


0.47 ±0.20 


0.42 


0.42 


JH 


IIP 


+30 


23.3 


22.5 


0.27 ±0.03 


0.36 


0.22 



a Relative to the epoch of B-band maximum 

b Redshift of the host galaxy, as determined by the Sloan Digital Sky Survey (SDSS) pipeline for 
SDSS galaxies 

c SNLS calculated redshift using broadband colors of SN 
d SNID redshift from our LDSS-3 spectra 
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Fig. 1. — Rost-frame wavelength spectra of supernova types, as well as an AGN template: SN la template (Hsiao et al. 2008), SN 1981b 
(Branch et al. 1983), SN 1999em (Baron ct al. 2000), SN 19941 (Clocchiatti ct al. 1996) and AGN LINER template (Francis et al. 1991). 
The spectra arc calibrated on the y-axis by 2.5 X log(Fjv/i 7 B), where Fjv is the flux in a one-angstrom narrow-band filter and Fg is the 
flux in a 500 A filter centered around the narrow filter. Offsets, shown by the solid lines, are added to the values. Compared to spectra 
like that of the SN IIP, or the AGN, the SNe la has very deep and wide features. 
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-1.0 L , , i , , , i , , , i i , I L i i i i i i i i i i i i i 

-0.2 0.0 0.2 -0.2 0.0 0.2 

^"^temp 

Fig. 2. — (left) The correlation function of a high SNR spectrum with a SN la template (solid) and SN II template (dotted). The peak 
value of the correlation curve allows one to determine the true redshift, in this case z — ztemp = 0. (right) The correlation function of 
a low SNR spectrum with a SN la template (solid) and SN II template (dotted). Here, a discrete set of redshifts, z — ztemp = 0, and 
z — ztemp = ±0.2 is preferred. The marks in the right panel demonstrate how with two measurements of the correlation curve, one can 
discriminate the type of SN, and for SNe la, determine its redshift to <r z S3 0.01. 
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Fig. 3. — The fourier transform of a SN la (z=0, peak brightness) spectral template (solid-triangles) and SN II (z=0, peak brightness) 
spectral template (dashed-squares) in the top panel. The SN la spectrum shows the majority of its correlation amplitude at k = 6, or 
~ 500 A, while the correlation amplitude for the SN II spectrum at k = 6 is much weaker. In the lower panel, a sine curve (dashed line) 
with wavelength 500 A is overplotted on a SN la (solid) and SN II spectra (dotted). 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



z=0.3 




Fig. 4. — Fourier transforms of SNe la spectra in the observed frame at redshifts 0.3 < z < 0.7. For each redshifted spectrum, there is 
a high correlation amplitude at k = 4, which is equivalent to a wavelength separation of 775 A. The peak wavenumber depends on the 
redshift z a such that k p = 6/(1 + z s ), but changes little over a modest fractional change in (1 + z a ). 
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Fig. 5. — The optimized three-tooth filter design, with one filter portrayed as a solid line, and the other as a dashed line. The g' and r' 
Sloan filters are seen in gray. SN la spectra at z = 0.38, 0.58 are overplotted to illustrate how the SNACC filter passbands correlate with 
the SNe la spectral features. A SN IIP spectrum at z = 0.50 is also shown to demonstrate little correlation between the SN II features and 
the filter teeth. 
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Fig. 6. — The CC — (g + r) diagram for the prototype design of the SNe la (outside) and the SNe II (inside). The SNe II redshifts cluster 
in the center of the diagram, and the SN la redshifts have an average radius m 0.25 mag. For measurement errors <ri 2 = 0.04, the SNe la 
are separated from the SNe II in CC — (g + r) space by 3 — 5 X (Ti,2, allowing for a SNe la sample purity of ~ 98%. 
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Fig. 7. — Histogram of redshift estimation with only a photometric redshift (dotted), only SNACC filters (solid) and the combination of 
SNACC filters with a photometric redshift (dashed). The tails in the population (3% for SNACC+Photometric, 45% for just SNACC) of 
redshift errors at ±0.2 are due to the degeneracy of the redshift in the correlation curves. Overall, the combination of SNACC filters with 
the photometric redshift yields a main redshift error of o z ~ 0.01. 
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Fig. 8. — Effects of reddening, stretch and age on the resulting fraction, Yj a , of SNe la correctly classified (type) for various redshifts. 
Reddening has very little effect on classification, except for the lowest redshifts, while stretch has a similarly small effect except for the 
highest redshifts. The filters were optimized for use around peak (±5 days), and the classification is generally > 90% this range. 
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Fig. 9. — Effects of ignorance of reddening, stretch and age on the errors in determined redshift. Variations in reddening and stretch 
have a very small effect, producing redshift errors < 0.02. The filters were optimized for use around peak (±5 days), and the redshift error 
for this age range is mostly < 0.02. Knowledge of the reddening, stretch and age would eliminate these errors. 
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Fig. 10. — The bandpasses, Iridian's expected design, and Iridian's actual design, as well as the strongest skylines ([OI], Na D with 
relative strengths as found at Magellan) in the optical range. Iridian was able to build filters with bandpasses of an accuracy better than 
5% of the specified widths, heights and positions of the teeth. 
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Fig. 11. — The calibration of the stars (crosses) in the G8V field with the set of Gunn-Stryker stars (diamonds), after the CC — (g + r) 
offset is incorporated. The correct CC — (g + r) offset is found by fitting all of stars in the field with g — r < 1.0 to the line formed of the 
Gunn-Stryker stars, and the calibration error is due to the differences between the set of stars in the field and the Gunn-Stryker set. For 
this set of stars, the CCi — (g + r) offset is 0.33 ± 0.01 mag and the CC2 — (g + r) offset is 0.24 ± 0.005 mag, where the errors here are 
folded into the entire photometric error of the classification process. Once the stars have been calibrated, the CC — (g + r) value of the 
G8V star can be placed in the CC — (g + r) diagnostic space. 
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Fig. 12. — Similar to Fig. |12| the calibration of the stars (crosses) in the M4V field with the set of Gunn-Stryker stars (diamonds), after 
the CC — (g + r) offset is incorporated. For this set of stars, the CC\ — (g + r) offset is 0.33 ± 0.015 mag and the CC2 — {g + r) offset is 
0.22 ± 0.01 mag. 
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Fig. 13. — Combining Fig. |11| the magnitudes of the stars (crosses) in the G8V field in CC — (g + r) space, along with the Gunn-Styker 
stars (diamonds) and SNe la (triangles). The photometric value (with error bars) of the G8V star is within one standard deviation of its 
synthetic position (highlighted square). 
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Fig. 14. — Combining Fig. |12| the magnitudes of the stars (crosses) in the M4V field in CC — (g + r) space, along with the Gunn-Styker 
stars (diamonds) and SNe la (triangles). The photometric value (with error bars) of the M4V star is within one standard deviation of its 
synthetic position (highlighted square). 
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Fig. 15. — Observed Spectra for NU, OE and JH with SN templates over-plotted. Using the SNID program, the supernovae are classified 
as: NU, a SN la with z = 0.56 ± 0.01; OE, a SN la with z = 0.42 ± 0.01; JH, a SN IIP with z = 0.22 ± 0.01. These correlations have SNID 
rlap values of 0.63, 0.68 and 0.74 respectively. 
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Fig. 16. — NU, OE and JH in CC — B space (with error bars), as well as their synthetic positions (highlighted symbols). The synthetic 
positions of the SNe are within the standard deviations of the photometric positions. The distribution of SNe la in CC — B space appears 
different than from Fig. [6] because here the Megacam sloan g' and r' filters make up the broadband filters. From the positions of the SNe 
in CC — (g + r) space, one can find the type of SN, and if la, its redshift. 
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Fig. 17. — The normalized probability distribution for both redshift and type (SN la, SN II) for NU, OE and JH from the chi-squarcd 
diagnostic (Eq. 8). The spectroscopic redshifts are also marked with a solid vertical line at its redshift value. For NU, the probability 
distribution for SNe la at different redshifts is a gaussian centered at z = 0.56 with a z = 0.01. The peak is ~ 20x greater than the 
probability (< 0.01) of the fiat SNe II distribution. For OE, the SN is clearly la, but there are degeneracies in the probability distribution 
due to the large photometric error and OE's position in CC — B space. For JH, there is no likely SN la candidate. 
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Fig. 18. — The optimized four-tooth filter design, with one filter portrayed as a solid line, and the other as a dashed line. The g' and r' 
SuprimeCam filters are seen in gray. SN la spectra at z = 0.38, 0.58 are overplotted to illustrate how the SNACC filter passbands correlate 
with the SNe fa spectral features. A SN IIP spectrum at z = 0.50 is also shown to demonstrate little correlation between the SN If features 
and the filter teeth. These four-tooth filters have an equivalent width of ~ 600 A. 
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Fig. 19. — The CC — B diagram for the four-tooth filter design of the SNe la (outside) and the SNe II (inside). While the SNe la radius 
is slightly smaller (0.20 mag), the cluster of SNe II is tighter around the center, and there is a similar separation between the SNe la and 
SNe II as from the three-tooth case, except here there is no local degeneracy of the SNe la. The SNe la purity and redshift determination 
for the three-tooth and four-tooth cases are approximately equal. 



